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Abstract 

Many models of Beyond the Standard Model physics involve heavy colored fermions. We study 
models where the new fermions have vector interactions and examine the connection between 
electroweak precision measurements and Higgs production. In particular, for parameters which 
are allowed by precision measurements, we show that the gluon fusion Higgs cross section and the 
Higgs decay branching ratios must be close to those predicted by the Standard Model. The models 
we discuss thus represent scenarios with new physics which will be extremely difficult to distinguish 
from the minimal Standard Model. We pay particular attention to the decoupling properties of 
the vector fermions. 
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I. INTRODUCTION 



The Standard Model of particle physics has a remarkable body of experimental support, 
but the Higgs boson remains a missing ingredient. Precision electroweak measurements 
suggest that a Standard Model Higgs boson must be lighter than ~ 145 GeV and 
recent measurements from the LHC exclude a Standard Model Higgs boson in the range 
129 GeV < Mu < 600 GeV |3j. Preliminary measurements suggest a light Higgs boson in 
the mass region Mh ~ 125 GeV {3, 4]. Should this putative Higgs signal be confirmed, the 
pressing issue will be understanding its properties. 

For all Higgs masses, gluon fusion is the dominant production mechanism at hadron 
colliders and the production rate is well understood up to NNLO in QCD {5,6]. Theoretical 
uncertainties from renormalization/factorization scale choices and from the choice of parton 
distribution functions are also well understood \l- 10|. The total rate, however, is sensitive 
to the existence of colored particles which couple to the Higgs boson. Beyond the Standard 
Model physics can potentially have a large effect on the Higgs boson production rate, making 
this a window to high scale physics js[ [ill 13|. 

The effects on Higgs production of squarks, Kaluza Klein colored fermions, color octet 
scalars, fermionic top quark partners and 4 th generation fermions (among many others) have 
been extensively studied. The simplest possibility for new heavy fermions is to form a chiral 
heavy new generation which, except for masses, is an exact copy of the known generations. 
After careful tuning, it is possible to find combinations of 4 th generation fermion masses 

nn 

which are permitted by precision electroweak measurements [14J, H5[ and are not excluded by 
direct searches. Since a chiral 4 th generation quark is assumed to couple to the Higgs boson 
with a strength proportional to its mass, heavy quarks do not decouple from the production 
of the Higgs boson (and in fact increase the rate by a factor of ~ 9). the existence of a 
4 th generation of fermions would exclude a Higgs boson mass up to Mu ~ 600 GeV [l6] 
regardless of the fermion masses. 

In this paper, we study the effect of heavy vector quarks on Higgs boson production 
and study both the case of an isospin singlet top partner and an isospin doublet of heavy 



fermions. A vector singlet top partner arises naturally in Little Higgs models [17H23| , where 
the couplings to the Higgs boson of the top quark and its fermion partner are fixed in 
such a manner as to cancel their quadratically divergent contributions to the Higgs mass 
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renormalization. Top-quark fermion partners are also found in top color 24j, |25) and top 



condensate 



models where there is a natural hierarchy of scales such that the top 



partner obtains a large Dirac mass. Light vector fermions instead typically appear in com- 



posite Higgs models 



3J33 



33|. Our results are general enough to be applied to any of these 
models and hence represent a simplification of results which have previously been presented 
in the context of very specific scenarios. 



A study of the S, T, U parameters and the Z — > bb decay rate |34j . [35| restricts the 
allowed parameter space for heavy vector fermions. However, vector fermions have interest- 
ing decoupling properties as the mixing with the Standard Model fermions becomes small, 
which makes a large region of parameter space experimentally viable. Vector fermions which 
couple to the Standard Model fermions and Higgs boson can be SU(2)l singlets with the 

same hypercharge as the Standard Model right-handed quarks, doublets with Y = Ysm = h 
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or Y = Ysm ± 1, or triplets with Y = Y SM ± | [36j. We consider the "Standard Model- 
like" case with either a heavy fermion singlet of charge 2/3 or a doublet with the Standard 
Model assignments of hypercharge. We compute the NNLO prediction for Higgs production 
for the allowed parameter region of these models and quantify the allowed deviation from 
the Standard Model prediction. The new features of our study include up-to-date fits to 
precision electroweak measurements in models with vector fermions, and an analysis of the 
the resulting consequences for Higgs boson production at NNLO in perturbative QCD. 

II. THE MODELS 

We consider models with additional vector-like charge 2/3 quarks, T Q , and charge -1/3 
quarks, B a , which mix with the Standard Model-like third generation quarks. For simplicity 
we make the following assumptions: 

• the electroweak gauge group is the standard SU{2) L x U(l)y group; 

• there is only a single Standard Model Higgs SU(2)l doublet, 




with 0° = 
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• we neglect generalized CKM mixing and only allow mixing between the Standard 
Model-like third generation quarks and at most one new charge 2/3 quark singlet or 
one new SU(2)l quark doublet. We do not consider fermions in more exotic represen- 
tations. 

The Standard Model-like chiral fermions are 

^=[ 72 |, 72, Bi, (2) 




with the Lagrangian describing the fermion masses 

- C S M M = \^\EB\ + + h.c. , (3) 

and H = ia 2 H*. 

The models we consider are: 





singlet fermion model: add a vector SU{2) L quark singlet of charge 2/3, 7^ and Tj 
doublet fermion model: add a vector SU(2)l doublet of hypercharge 1/6, 
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III. EXPERIMENTAL LIMITS ON TOP PARTNER MODELS 
A. Limits from and A\, 

Data from LEP and SLD place stringent restrictions on the couplings of the fermionic 
top partners. The top partners mix with the Standard Model-like top quark and contribute 
at one loop to processes involving bottom quarks, especially Z — >■ bb and Af,. The neutral 
current couplings to the bottom can be parametrized by the effective Lagrangian 



C NC = ^-Zj-f 
cw 



9 { + ^)(^Wj£ + ^ A/1 



/, (5) 



where the Standard Model couplings are normalized such that g* L = T( — QfS^r, 
g f R = -Q f s 2 w , with s 2 w = sm 2 6 w = (e/g) 2 = 0.231 Q and T* = ±1/2. The couplings 
5g^ L R = Sg^^ 1 + 5g-[ R contain both the Standard Model radiative corrections, Sg^ 3 ^ 1 , and 



95% Confidence Level Limits 



-Rb 




FIG. 1: Allowed 95% confidence level regions from the simultaneous fit to Rb and Ab (red shaded), 
Rb alone (between solid black lines), and Ab alone (between dashed blue lines). 



the new physics contributions, 5gt R . The Standard Model contribution from top quark 



loops is well known 38M40| . and in the limit m t » Mz it is given by 



b.SM G F ml 



(6) 



v/2 87r 2 

The dominant effect of new physics in the b sector can be found by assuming that 8g b L 
and 5g b R are small and approximating 4l|, |42] 

V{Z -> 66) 



Rb 



Rf M l 1-3.57^ + 0.655^ 



T(Z —7- hadrons) 



Af M U-0.315g b L - 1.728g 



R ( 7 



(7) 



(5g b L r + «) 2 

where R% M and Af M are the theory predictions including all radiative corrections. The 



positive contribution to 8g b £ SM from the top quark has the effect of reducing both Rf M and 

M 

The 95% confidence level ellipse for new Zbb couplings is shown in Fig. [Hand is obtained 
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using the Particle Data Group results [37|, |43 ] 



R e b xp = 0.21629 ±0.00066 



R: 



SM 



0.21578 ±0.00005 



A e b xp = 0.923 ±0.020 
Af M = 0.9348 ± 0.0001 



(8) 



If 5g R = 0, the 95% confidence level limit from the fit to Ab and Rb is 



- 0.0027 < 5q b r < 0.0014. 



(9) 



Similarly, if 5g b L = 0, the 95% confidence level limit from the fit to Ab and Rb is 



- 0.0066 < Sg b R < 0.0148 . 



(10) 



B. Limits from the Oblique Parameters S, T and U 



The new quarks contribute at loop level to the vacuum polarizations of the elec- 
troweak gauge bosons U^ Y {p 2 ) — ^xy{p 2 )9^ u + B X y(p 2 )p^P u , with XY = 77, 7Z, ZZ and 



W + W~ 
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and Takeuchi 



45 [.These effects can be parametrized using the S, T and U functions of Peskin 





^s w c w J 




M z \ 














aUp 





n zz (Mf)-n zz (o)-n 77 (Ml) 
n^z(o) 



' w 



c w s w 



Ml 



M 2 W 



-w 



Il zz (M 2 z )-U zz (0) 



Ml 



n 7Z (Mf) 2 n 77 (Mf) 

-ZSwCw ttt^ S T . 



Ml 



■>w~ 



Ml 



(11) 



Any definition of sw can be used in Eq. [TT] since the scheme dependence enters at hi ghe r 



order. Since these parameters are well-constrained by LEP and LEP2 measurements 
they set stringent limits on the masses and couplings of the new quarks. 
We use the fit to the electroweak precision data given in Refs. 2j, |47I |. 

AS = S-S S m = 0.02 ±0.11 
AT = T -T SM = 0.05 ±0.12 
AU = U-U SM = 0.07 ±0.12, 
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46|, 



(12) 



with reference Higgs and top-quark masses Mn,mi = 120 GeV and m tjre f = 173.1 GeV. The 
associated correlation matrix is 



Pij 



( 1.0 0.879 -0.469 ^ 
0.879 1.0 -0.716 
-0.469 -0.716 1.0 



The Ax 2 is defined as 



A X 2 = ^(AX 4 - AX l )(a 2 ) i r. 1 (AX J - AX 



31 i 



(13) 



1,3 



where AX, = AS, AT and AU are the central values of the fit in Eq. [TJ, AX = X -Xf M = 
ASf, ATp and AUf are the contributions to the oblique parameters from the new fermions 
and crfj = UiPijUj , <Ji being the errors given in Eq. [12j A 95% confidence level limit in a 
three-parameter fit corresponds to A% 2 = 7.82. 

Since we consider primarily M H = 125 GeV, we need to add the Higgs contributions 1 

1 



AS H 
AT h 
AU h 



12tt 



log 



M 2 H 



-3 



' Ml 



log 



M 2 H \ / 'M 2 7 



Ml rei 



(14) 



C. Other Experimental Limits on Top Partner Models 



Both ATLAS 50|, |5l| and CMS 52|, |53j have searched for direct pair production of new 
heavy fermions. For charge 2/3 top-like quarks decaying with 100% branching ratio to Wb, 
CMS excludes masses below 557 GeV at 95% confidence level, while ATLAS sets an upper 
bound of 404 GeV. CMS dedicates a specific analysis to pair-produced vector quarks of 
charge 2/3 decaying entirely to Zt, excluding masses below 475 GeV js4|. For charge -1/3 
quarks, assuming 100% branching ratio to Zb, ATLAS excludes masses below 358 GeV for a 
vector singlet, while CMS excludes charge -1/3 quarks decaying with 100% branching ratio 
to Wt below 611 GeV. These limits are not directly applicable to our models, since the 



1 Our fits include the exact results for the Higgs contributions, which can be found in many places including 
Ref. 



48( and the Appendix of Ref. [4' 
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branching ratios of the new heavy fermions to Standard Model particles are degraded by 



mixing angles and the limits therefore weakened [36ll55ro9j|. Our results on Higgs production 
are rather insensitive to the masses of the new top partners and we typically assume masses 
of the TeV scale. 

In principle, there are also limits on heavy charged fermions which mix with the Standard 
Model third generation quarks coming from K, B and D rare processes. For TeV-scale 
masses of the new fermions and small mixing parameters (which we will see in the next 
section are required by limits from oblique parameters, Rb and Ab), the constraints from 



rare processes are not restrictive 



57|, 



59 



461]. 



IV. SINGLET TOP PARTNER MODEL 



Little Higgs models [17H23|. topcolor models [24j, |25j and top condensate models (26M28 ] 
all contain a charge 2/3 partner of the top quark, which we denote by T 2 . We consider a 



general case with a vector SU 
Model-like top quark 



19 



36 



2)j, singlet fermion which is allowed to mix with the Standard 
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63J. The mass eigenstates are b = B 1 , t and T, where 



b and t are the observed bottom and top quarks. Thorou. 



measured mass values m& = 4.19 GeV, m t = 173.1 GeV 
charge 2/3 can be found through the rotations 

tL,R 

Tl,r 
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i this paper we will use the 



64|. The mass eigenstates of 



T 1 

1 L,R 



L,R 



T 2 



(15) 



The matrices Ut R are unitary and r = ^r 1 ^- The mixing matrices are parametrized as 



L,R "—""-j ^ - L,,n — 2 

cos 6l — sin#L 
sin^L cos#l 
We abbreviate cl = cos 8l, sl = sinflj,. 
The most general fermion mass terms are 

-i 



J J 1 

U R 



cos 9 ft — sin 9r 
sin 9r cos 9r 



(16) 



-C 



-c s m m + wiifrz + ktIu + + h.c. 



M,l 



Xl 



UliM^ + hH^jUx 



f.f 



X t R + X 1 V -^B 1 L B R + h.c. 



where 



t 

(i) 



\ V \ V 

A2 75 A3 7I 



H m = — 




(17) 



(18) 



The resulting mass eigenstates are 

M t,diag = | " I _ ^ 




One can always rotate T 2 such that A4 = 0. Since A4 can be rotated away, the model has 
four free parameters. Alternatively, it is always possible to rotate T R such that sin6 R = 0, 
because only the Standard Model-like left-handed doublet mixes to the singlet with a 
Yukawa term 2 . Therefore the couplings only depend on 9l, which we will take as one of the 
four physical parameters along with mj, (physical mass of the charge -1/3 quark), m t and 
M T (physical masses of the charge 2/3 quarks). 

The physical masses and mixing angles are found using bi-unitary transformations, 

= UftMftMfautf. (20) 

It is straightforward to find the mass eigenstates and mixing angles, 

2A4A5 + v 2 A 2 A 3 



tan(2#R) 
tan(20 L ) 



v 

m t M T = — I A 2 A 5 - A3A4 



A 2 - A| + y (A§ - A|) 
V2v(X 2 X 4 + A3A5) 

A 2 + A| - \ + A|) 

v 

7! 



M| + m 2 = y(A 2 + A 2 ) + A 2 + A 2 . (21) 
From Eq. [TTJ the couplings to the Higgs boson are 

CI = c tt t L t R h c TT T L T R h c tT t L T R h c Tt T L t R h + h.c. , 22 

V V V V 

where 

ctt = c\ , c T t = s\ , c fT = c Tt = s L c L . (23) 
These relations can be easily derived by noticing that 

vH\ l)M = M[ 1)jks 8 kl , (24) 



2 We do not perform any of these rotations here, and the formulas in this section hold for the arbitrary 
Yukawa couplings of Eq. [T7] 



9 



yielding for the physical Higgs couplings 



Hi),* 
V 1U lJki 



R,rs 



TjtA 

u R,sj 



where the index k is not summed over. 

The charged and neutral current interactions are 



cc 



and 



NC 



C W 



i=t,T ^ 



=1,2 



u 



I ) 7 A + h.C. 



(9l + Sgl)P L + (g R + 5~g R )P R 



fi 



/i 



3 ( ' 



where R = 5g l l s R + <5(?£ R contains both the Standard Model contribution from top 
quark loops and the new physics contributions. The new physics couplings arising from the 
interactions of the top partner singlet are 



(25) 



(26) 



(27) 



2 ' 

,y - o 



i,j =t,T. 



HI = ^ 



(28) 



Finally, the unitarity bound from the scattering TT — > TT is modified from the Standard 



Model limit and becomes 



65 



66) 



(Unitarity Bound) 



< 



550 GeV 



(29) 



This class of models therefore allows quite heavy T quarks without violating perturbative 
unitarity. 



A. Experimental Restrictions on Singlet Top Partner Model 



Using the results given above, it is straightforward to compute the contributions to the 
S,T and U parameters in the singlet top partner model. Expressions for the gauge boson 
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two-point functions for arbitrary fermion couplings are given in the appendix 2l|, [67J, |68[ , and 
we assume Mt » Mz- Subtracting the Standard Model t — b loops, the new contributions 
are 



ATp 
AS F 



rp 2 

J-SM S L 



:i + cj) + sjr + 2ci 



r — 1 



log(r) + O 



Ml Mi 



m 2 ' M T ' m 2 



N, 



C „2 



18tt 



sU\og(r)+c L 



5(r 2 + l)-22r 3(r + l)(r 2 - 4r + 1) 



(r - l) 1 



+ 



[l-rf 



log(r) 



( 



Ml M| 



m 



\ m 2 ' M|, ' m 2 



AS F + A£/f 
where 



N, 



C „2 



9vr 



log(r) + O 



Ml Ml 



mt 



ml ' M T ' m 2 



M 2 
2" 



c 



and T. 



c 



SM 



16ns 2 w M 2 , 



(30) 



(31) 



is the t — b contribution to the T parameter in the limit of a massless bottom quark. We use 
M w = 80.4 GeV, Mv = 91.2 GeV [37]. Eq. M agrees with the m b ->■ and M z « M T ,m t 



limits of Refs. 



69 



7l). 3 For a top parter much heavier than the top quark, Mt » m t , the 



oblique parameters take simple forms, 



ATp(approx) 
ASV(approx) 



rp 2 

J-SM S L 



N, 



18tt 



C s 2 



-(l + c|) + slr + 2c 2 log(r) 
5 c 2 +(l-3c 2 )log(r) 



(32) 



One would expect decoupling to occur for a very heavy vector top partner, i.e. for 
r — > 00. From Eqs. [301 E21 instead, the effects of the top partner on the oblique parameters 
vanish only in the limit sl — > 0. This can be understood inspecting the mass matrix flTSJ): 
to obtain decoupling both the Yukawa interactions and the off-diagonal terms need to be 
small, X2V, X^v, A4 << A5. In this limit 

v\ 3 



A2A4V 



V2M 7 



(33) 



V2M T v-.w :r 

and the top partner effects vanish for large Mt, as expected. 

In Fig. [2] we show the oblique parameters for a fixed Mt = 1 TeV. It is clear that the 
approximate forms of Eq. [32] are excellent approximations to the complete results for mass 



3 There is a typo in Eq. 88 of Ref. [ic]] , where there is a 2 instead of the factor 22 in Eq. [30l 
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Singlet Fermion Model 



\ M = 1 TeV - 
^ — 

: - 




i 
~- 

,J 

: 








- AT F 


- 




-- AT F (approx) 




-1 r 

/ /y 

/ A 

Iff 


— as c 

F 






-- AS F (approx) 






— all 

F 






-- AU F (approx) 











0.1 0.2 0.3 0.4 

sin 8 L 



FIG. 2: Contributions to ATp , ASf, AUf from a singlet top partner as a function of sin#p 
for fixed My = 1 TeV. The results of Eq. [32] in the limit Mt » nit are shown as ATp(approx), 
A,Sp(approx) and At/ir(approx). 



Singlet Fermion Model 



= ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 
sin 6 L = . 1 


, i | i i i | 

— AT p 

. .. AT F (approx) 

— AS C 

. .. AS F (approx) 

— ALL 

F 

■ ■■■ AUp(approx) 









800 1000 1200 1400 1600 1800 2000 

M T (GeV) 

FIG. 3: Fermion contributions to ATp , ASf, AUf in the singlet top partner model for 
fixed sm6i = 0.1. The dotted lines represent the approximate results from Eq. [32] in the limit 

M T » m t . 
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Singlet Fermion Model 



sin 6 L = .05 
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FIG. 4: Same as Fig. [3] for a smaller sin 9l = 0.05. 



2000 



values of this order. The largest contribution is to ATp due to the large isospin violation for 
non-zero sin^p. In this case, the isospin violation is manifest in the result that AUp > ASp. 
The new physics effects vanish as sl — > and we recover the Standard Model couplings. 
The oblique parameters for fixed sin^ are shown in Figs. [3] and H] as a function of Mt- As 
we argued, the decoupling does not occur for large M T , but requires sl —> 0. 

The parameter space allowed by a fit to the oblique parameters can be found using the 
results of Section IIHI Fig. [5] shows the 95% confidence level upper bound on the mixing 
angle sin 6l as a function of Mt for a light Higgs boson. For a heavier Higgs boson, it is 
possible to use the positive contribution to T from the top partner to compensate for the 
negative contribution from the heavy Higgs, as shown in Fig. El In this case, a minimum 
degree of mixing is required, since such a heavy Higgs boson is excluded by the electroweak 
fit in the three-generation Standard Model. The heavier the Higgs boson, the smaller the 



range of sin^p allowed. 
(M H ~ 800 GeV) in Ref. 



his situation was explored for an extremely heavy Higgs boson 



72]. 



The mixing in the charge 2/3 sector also affects Rb- In the limit m& — > and neglecting the 
small correlations between Rb and the oblique parameters, only 8g\ is affected by the singlet 
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95% Confidence Level Upper Bound from STU 




FIG. 5: 95% confidence level upper bound on the mixing angle sin^L in the singlet top partner 
model from experimental restrictions on the S, T and U parameters. 




FIG. 6: 95% confidence level bands allowed by a fit to S, T and U for the mixing angle sin#£ in 
the singlet top partner model. The regions allowed are between the curves corresponding to each 
value of Mh ■ 
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top partner. Its contribution to 5g b L can be found from the general analysis of Ref. [7c 



9 2 2 



647T 2 



s L [ fifaa/) + c L f 2 (x,x') 



(34) 



where x = mf/M^ and x' = M^/M^. The Standard Model top contribution has been 
subtracted following the definition of Eq. [5j In the heavy mass limit, x,x' » 1, 



fi{x,x') 
f 2 (x,x') 



x 



/ X 

x + 3 log — 

\ X 

2xx' 



-x 



X + 



X' 



X 



log — 



X 



X 



(35) 



The contribution to bg\ from top singlet partners is shown in Fig. [7] as a function of for 
fixed sin^L, along with the 95% confidence level region allowed from the fit of Section l"HI Al 
We use the exact one-loop calculation of 5g b L in the top singlet partner model, following 
Refs. 0,13 • 

The relatively large contributions from R b can be understood by looking at 
the leading terms for m t , Mt » Mw, 



r b x b,SM 2 

09l = °9l s 



:i + c 2 L ) + s 2 L r + 2cl- 



1 



log(r) 



(36) 



Again we see that the heavy top partner decouples only when the parameters in the mass 
matrix are such that sl ~ Furthermore, its contributions to 8g b L and to the T parameter 
(Eq. I32i) are both positive and strongly correlated. A large, positive contribution to 5g b L from 
the singlet also implies a large contribution to the T parameter. This correlation was already 
pointed out in the context of composite Higgs models in Refs. 



31 



32 



75J. 



Combining the new physics contribution with the Standard Model top quark contribution, 



Sg b L = 5gr m + Sgl = 5g 



b,SM 



b,SM 



cl + slr + 2s\cl 



r — 1 



log(r) 



(37) 



and the net effect of the top partner is to increase Sg b L and hence reduce R b . 

A comparison of the maximum value of sin 9l allowed by the fit to R b and A b and by the 
experimental limits arising from the fit to S, T and U (with Mh = 125 GeV) is shown in 
Fig. El where it is apparent that the most stringent restrictions in the top partner singlet 
model come from the oblique parameters. 
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Singlet Fermion Model 

95% Confidence Level Limit from Z— >bb and A, 





— sin 6 L = .05 
... sin8 L = 0.1 
-- sin8 L = 0.2 




model. The 95% confidence level allowed region is shaded in light blue. 



95% Confidence Level Upper Limits 



— Limits from STU with M =125 GeV 

ri 

- - Limits from R, 




M T (GeV) 

FIG. 8: Maximum allowed sinfjf, in the top partner singlet model from oblique measurements 
(black solid) and Rb (red dashed). 
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V. TOP PARTNER DOUBLET 



A. Model with Top Partner Doublet 



In this section, we consider a model which has in addition to the Standard Model fields 
a vector SU(2)l doublet [361. |76| . 



^R 



n 



(38) 



The most general fermion mass terms allowed in the Lagrangian are 



C 



M.2 



R 



'R 

U b L M b (2) U b J 



Xr + h.c. 



where Xl r are given by Eq. [15] and 



(39) 



Xl,r 




U 



L,R 



Kr, 



(40) 



We can always rotate "02 such that Ag = 0. So without loss of generality 



M\ 2) 



^ 
X7 72 As 



M, 



(2) 



^7=2 



Because of the SU(2) symmetry 



(41) 



(42) 



Diagonalizing the mass matrices now requires four unitary mixing matrices, U l L , U R , U R , 




(43) 



There are five input parameters in the Lagrangian. We will take the five independent 
physical parameters to be the physical masses, m t , Mt and m&, Mb (with m t and m;, being 
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the mass of the Standard Model-like fermions) and the right-handed mixing angle in the 
charge -1/3 sector, 9 R . It is straightforward to find relationships among the mixing angles: 



(sin9 R f = - (8in ^ )2 



^sin^) 2 + (cos^) 2 9r 

(sin^f = ^ I 

(sin#5,)2+(co 8 »?,)^ 

\2/„,2 T\/r2\ i A/f2 /„• nt \2/2, i\/r2\ , 71 ^2 



(sm9 b R Y(m 2 b - M B ) + M B = (sin 9 R ) {m t — M?) + Mf. . (44) 

For small mixing, the left-handed angles are always suppressed by the heavy mass scale 
relative to the right-handed angles of the same charge sector, 

mf 



■ 2 nt,b "H,b . 2 Q t,b / ,r\ 

sm l ~ TJT- sm e R ■ ( 45 ) 

T,B 



If the mass splitting between B and T, 5 = My — Mb, is small, 4r- << 1, 



The charged current interactions are 



+ ^r1^ iPr > + h.c 



where 



v/2' 

-^^jxlT^i^xi + + h.c. , (47) 

a-=( 01 ]. (48) 

The neutral current interactions are given by Eq. [27J with 

Sgl = 6g% = 0, i,j = t,T,b,B 

5g R = Tl sin 2 9% , 5g T R = T 3 T cos 2 9 R , dg^ = - - sin ^ cos 9 R 

6g b R = Tl sin 2 9 b R , ^ = T 3 b cos 2 ^ , <fo*? = \ sin ^ cos 9 b R . (49) 

In the doublet top-partner model, only the right-handed Standard Model-like singlets 7^, 
Bf have Yukawa type mixing with the new quarks, and therefore only the interactions in 
the right-handed sector are modified. Finally, the Higgs couplings are given by 

A = -°tt — t L t R h - c TT T L T R h - c tT — t L T R h - c Tt T L t R h + 

V V V V 

-c bb —b L b R h - c B B — -B L B R h - c bB —bLB R h - c B b — -B L b R h + h.c. , (50) 

v v v v 
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where 



c a = cos 2 9 R , c TT = sin 2 9 R , c tT = c Tt = sin 9% cos 9% 



c bb = cos 2 9 b R , ess = sin 2 9 b R , c bB = c B b = sin 9 b R cos 9 b R . (51) 
The derivation of these couplings follows the lines of Eq. just with 

vH\ 2)M = M[ 1)M 5 sl . (52) 

B. Experimental Limits on Doublet Fermion Model 

The decay Z — > bb puts a strong restriction on sin 9^ since mixing in the right-handed 
6-quark sector contributes to 5g b R at tree level, 

<^ = -~(sin^) 2 - (53) 
From Eq. [101 the mixing angle in the right-handed b sector is highly constrained, 

| sin < 0.115. (54) 

The contribution to 6g b L in the left-handed sector occurs at one-loop. Subtracting out the 
Standard Model contribution, in the limit x, x' » 1, the approximate result in the doublet 



fermion model is 



8& = ^{sin 2 (0l - 9 b L )f 1 (x,x') + (sin^) 2 cos 2 (0i - 9 b L )f 3 (x, a;')} , (55) 



where f\(x,x ; ) is defined in Eq. [351 and 



tj / iAj \ I Jb | Jb \ iA/ , / tAJ 



A(«, *<) = -x + - ^1 + - J + - 3) ^ log^- j . (56) 

In Fig. [9l we scan over sin (9^ and and use the relationships of Eq. H31t o find the remaining 
parameters. We use the exact one-loop result for 5g^ following Refs. Q, and determine 
the 95% confidence level upper bound on Because of the tree level mixing in the b 
sector in this model, along with the relationships of Eq. 0H the heavy fermions are required 
to be approximately degenerate, as is clear from Fig. [9j This result is relatively insensitive 
to M T . 

Since sin 9 R is constrained to be quite small, we will consider the oblique parameters in 
the limit 9 R = 0. From Eq. 0U 9 R = implies 9\ = and the free parameters are m t ,mi,, 
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95% Confidence Level Limits from Z— >bb and A, 

b 

Doublet Fermion Model 




sin8„ 

K 

FIG. 9: Maximum value of allowed at the 95% confidence level from Rb and Af, as a function 
of sin 9%, in the doublet fermion model. 
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FIG. 10: Oblique parameters in the doublet fermion model in the limit @r = and mj — > as a 
function of M T . We fix 6 = M T - M B = 1 GeV. 
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FIG. 11: Maximum mass splitting 5 = Mt — Mb as a function of Mt allowed at the 95% confidence 
level from the oblique parameters in the doublet fermion model for 9 h R = . 



Mt and Mb- Our results will always be expressed in terms of 8 
9 R = 0, ui}, — > limit, the oblique parameters are well approximated by 69j, |76] 



M B . In the 



and 



ATp(approx) 
ASV(approx) 
AL r i?(approx) 

&9b = 5 9 



4T. 



SM 

Nc_^_ 

97T M T 

NcJ_ 
9vr M t 

5 



2 log(r) - 3 + 



5 log(r) — 3 



4 log(r) - 7 + 



41og(r) + 7' 



3 + 



6 log(r) — 17 



L,SM 

b M 



log(r) -4 + 3 



log(r) — 2 



(57) 



(58) 



It is apparent that decoupling of the effects of the new fermions occurs in the isospin conserv- 
ing limit, — » 0. Note also that 5gj^ changes sign for Mt — 7m f . The oblique parameters 
are shown in Fig. [10] for 5=1 GeV. As in the singlet model, ATp >> ASp, AUp. 

The limits coming from the oblique parameters are found from a global fit to S, T and 
U as described in Section HUB For 9 h R = 0, the 95% upper limit on the mass splitting, <5, is 
shown in Fig.[TT]as a function of My. For M# = 125 GeV and Mr ~ 1 TeV, the experimental 
constraints on the oblique parameters require 8 < 8 GeV. As shown in Fig. [TTJ it is possible 
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to compensate for the negative contribution to T from a heavier Higgs boson by a larger mass 
splitting 5, which generates a positive contribution to ATp. However, the limits from Aj, and 
Rt, in this model are much more stringent than those coming from the oblique parameters. 



VI. PHENOMENOLOGY 



by 



he new fermions affect the gluon fusion production rate, which at lowest order is given 



77 
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a(gg -»■ H) 



8 1- 



Ml 



(59) 



1024tto 2 

The sum is over t, b, T in the singlet fermion model and over t, b, T, B in the doublet model 
the Yukawa couplings normalized to the Standard Model values c qq are given in Eqs 
and l5ij and 



Fl/2 

f(r q ) 



AMI 
-Ml 



-2t„ 



1 + (1 - r q )f(r q 



sm 



log 



171 



Hr q >l 



if T q < 1 . 



(60) 



We compute the gluon fusion production cross section through NNLO using the program 
iHixs jsj]. iHixs allows the calculation of the cross section at NNLO for extensions of the 
Standard Model with an arbitrary number of heavy quarks having non-standard Yukawa 
interactions [13], and puts the predictions on a firm theoretical basis. At NNLO, there 
are contributions which mix quark loops of different flavors (e.g. t and T) and cannot be 
obtained by a simple rescaling of the lowest order cross section. We scan the parameter 
space allowed by the precision electroweak data as determined in the previous sections and 
discuss the maximum deviations from the Standard Model predictions. 

A. Top Partner Singlet Model 



The deviation from the Standard Model prediction of the NNLO Higgs production cross 
section as a function of the mixing angle in the top partner singlet model is shown in 
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FIG. 12: The ratio of the NNLO Higgs cross section in the top partner singlet model normalized 
to the Standard Model prediction as a function of s'mOi for Mh = 125 GeV and y/s = 8 TeV. The 
vertical line represents the maximum value of sin Ol allowed by electroweak precision measurements. 



Fig. [T2j The largest value of sin di allowed by the precision electroweak limits derived in 
Section IIV Al is also shown. As sin 6l increases, the mixing with the Standard Model-like 
top quark becomes significant, causing a suppression of the rate. This can be understood 
from the heavy mass limit (m t , >> of the lowest order cross section, where the 
gluon fusion rate scales as 

c7sm 60 m? L \ M*J 1 ; 

However, only the region to the left of the dot-dash line in Fig. [12] is allowed by the precision 
electroweak measurements, making the Higgs boson production rate in this model almost 
identical to the Standard Model rate. In contrast with composite [s^l or Little Higgs 59, 81] 
models, which typically have a sizeable reduction of the Higgs production cross section 
relative to the Standard Model, in models with vector fermions the suppression is negligible 
because of the decoupling properties discussed in the previous sections. The uncertainty on 
the Standard Model cross section coming from scale, PDF, and a s uncertainties is roughly 
15 — 20% so the extremely small deviation from the Standard Model prediction in the 
top partner singlet model is unobservable. The cross section for a heavier Higgs boson of 
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FIG. 13: The ratio of the NNLO Higgs cross section in the top partner singlet model normalized 
to the Standard Model prediction as a function of s'mOi for Mh = 300 GeV and y/s = 8 TeV. The 
region between the vertical lines represents the region of sin#£ allowed by electroweak precision 
measurements. 
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FIG. 14: Deviation of the H — >• 77 branching ratio in the top partner singlet model from the 
Standard Model prediction. 
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mass Mh = 300 GeV is shown in Fig. [131 In this case, there is a region of mixing angles, 
sinOi, which is allowed by the precision electroweak measurements. Again, there is a slight, 
but unobservable, suppression of the NNLO rate relative to the Standard Model rate. 

The loop-mediated Higgs decays to 77, Z7 and gg are also affected by the presence of 
top fermion partners. Fig. [TH shows the deviation of the branching ratio to 77 from the 
Standard Model prediction. For small mixing, this deviation is always less than one percent. 

B. Top Partner Doublet Model 

The deviation from the Standard Model prediction for the NNLO gluon fusion cross 
section for Higgs production in the top partner doublet model (computed using iHixs) is 
shown in Fig. [T5j Also in this case the maximum difference from the Standard Model in the 
allowed regions of parameter space (Fig. [9]) is always less than a few percent. This result 
can be understood by considering the heavy mass limit of the lowest order cross section for 

Doublet Fermion Model 

Vs=8 TeV, pp^H 

1 .02 1 1 1 1 1 1 1 1 1 1 1 1 1 




FIG. 15: The ratio of the NNLO Higgs cross section in the top partner doublet model normalized 
to the Standard Model prediction as a function of sm.9 b R for Mh = 125 GeV, yfa = 8 TeV, and 
M T = M B = 1 TeV. 
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the gluon fusion production of the Higgs, 



^Doublet / 1 i • 2 nb 

' 1 1 + Sill R 



&SM 



inh 7 M 2 H 2r- 1 2nh 

i + s n 2 ei f s n 2 e b R 

R 60 m 2 V r R 



26 25 r + 1 . 2 , 
sur 6>p 



(62) 



From the fits to Ab and the maximum value of sin 6^ is restricted to be 0.115, which 
implies 

^i£l < (i + sin 2 ^ 2 ~ 1.03 . (63) 

Similarly, the deviations from the Standard Model in the Higgs decays to 77, Z7 and gg 
and in H — > bb, which is affected at tree level, are not observable due to the small mixings 
and mass splitting allowed. 



VII. CONCLUSIONS 

We have considered the effects on the gluon fusion Higgs boson production at NNLO 
from heavy vector quarks of charge 2/3 and -1/3. Since the new quarks are vector-like, their 
couplings to the Higgs boson are suppressed by mixing angles relative to the Standard Model 
Yukawa couplings. These mixing angles are restricted to be small by precision electroweak 
measurements. The most stringent bounds come from the oblique parameters for a vector 
singlet top-partner, and from Af, and Rf, for an extension of the Standard Model with an 
additional vector doublet. Because of the small mixing angles allowed, in these models the 
Higgs boson production rate as well as its decay branching ratios are essentially those of the 
Standard Model. The scenarios we have presented will be extremely difficult to disentangle 
from the Standard Model without the observation of direct production of the heavy fermions. 
Vector doublet fermions with a non-standard hypercharge assignment are less restricted by 



precision electroweak measurements [36J] and the mixing angles between the t — b sector and 
the new fermion sector can be larger than in the cases we considered. However, even in 
this case, the low energy theorems for Higgs production require that the Higgs cross section 
approach the Standard Model result for large fermion masses. 

If a Higgs boson is found at the LHC, attention will turn to understanding its properties. 
By performing global fits to the measured rates, information can be gleaned from the various 
cross section times branching ratio channels. For a light Higgs boson, it is likely that the 
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dominant production channel will be gluon fusion, even in models with new physics. In 
this case, the rates are sensitive not only to a rescaling of the Standard Model couplings, 
but also to the effects of new particles which couple to the Higgs boson and contribute 
to the decay rates. Numerous preliminary attempts have been made to use current LHC 
data to discern differences from the Standard Model 
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901 ] . Our scenario with vector 



fermions demonstrates the difficulty of these indirect determinations of new physics - it is 
(un)fortunately not difficult to construct models which give Higgs signals indistinguishable 
from the Standard Model. 
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Appendix: Two Point Function for Arbitrary Fermion Coupling 



The contributions to the gauge boson two point functions from fermion loops 
parametrized by the interaction 



for V = W, Z, 7 are 

Nr. (2 
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67] 



£ = fi (c f L fp L + c f R fp R } 7/l / 2 v 



(64) 



n 



XY 



167T 2 1 3 V x Y RX Y 



m x + m\- — - ( A (mi) + A (m 2 ) 



2p 2 



2p 2 



+2m x m 2 (c f L fC s ^ + C f R fC f L ^ B, 



Q (mi,m 2 ,p 2 



(65) 



where 



and 



A (m) = 
/i(mi,m 2 ,p 2 ) = 



47r/i2V r(i + e) (i + iK 



+ 6 



- - fi(m u m 2 ,p 2 ) 
e 



dx log I x + 



m\(l — x) — p 2 x(l — x) 



(66) 



(67) 
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